GABOR REPRESENTATIONS OF EVOLUTION OPERATORS 
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Abstract. We perform a time-frequency analysis of Fourier multipliers and, more 
generally, pseudodifferential operators with symbols of Gevrey, analytic and ultra- 
analytic type. As an application we show that Gabor frames, which provide optimally 
sparse decompositions for Schrodinger-type propagators [14], surprisingly reveal to 
be an equally efficient tool for representing solutions to hyperbolic and parabolic-type 
differential equations with constant coefficients. In fact, the Gabor matrix represen- 
tation of the corresponding propagator displays super-exponential decay away from 
the diagonal. 



1. Introduction 

A Gabor system is obtained by fixing a function g G L 2 (K d ) and considering the 
time- frequency-shifts 

(1) 9m.n = e 27rinx g{x - m), (m, n) G dL d x (3Z d , 

with a > 0, j3 > 0. A Gabor system is a Gabor frame if there exist A, B > such that 

(2) A||/||i 2 <^|(/ !5mjn )| 2 <B||/||i 2 

m,n 

for every / G L 2 (R d ), see for example [H EHl EE]. With respect to fr ames of wavelets, 
curvelets and shear lets, in ([I]) dilations are then replaced by modulations. From the 
analytic point of view, this provides symmetric properties with respect to the Fourier 
transform and a simple treatment, cf. |21j . 

Gabor frames turned out to be appropriate tools for many problems in time- frequency 
analysis, with relevant applications to signal processing and related issues in Numerical 
Analysis, see for example [9], [34] . and references there. More recently, attention has 
bees addressed to the analysis of the partial differential equations [21 [31 01 [TUl [HI H21 
[131 EH E3 HO] , let us also mention the pioneering work [17] . 

Aim of this paper is to represent in terms of Gabor frames the solutions of the 
Cauchy problems for a class of evolution operators with constant coefficients, including 
hyperbolic and parabolic operators. We shall prove in particular a super-exponential 
decay away from the diagonal, for the corresponding Gabor matrix. 

To give a general setting for our results, let us denote by T the linear operator 
providing the solutions of a well-posed problem for a partial differential equation, in 
suitable function spaces. Generically, we expect that T is written in the form of 
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pseudodifferential operator (PSDO), or Fourier integral operator (FIO). Via its action 
on a fixed frame, in particular a Gabor frame, T can be seen as an infinite matrix, 
with entries 

(3) 2~m',n',m,n = {1 9m,ni 9m' ,n') ■ 

To deduce results of sparsity for this representation, one is then led to estimate the 
decay properties of T m i )n i mn , for large values of the indices m',n',m,n. 

In this order of ideas [5JE] obtained optimal sparsity for strictly hyperbolic problems 
in the frame of the curvelets. The operator T was a classical FIO with phase function 
homogeneous of degree 1 in the dual variables, cf. j27J Vol. 4], 

By parallel arguments the authors of |14j proved optimal sparsity for Schrodinger 
equations in the Gabor frames. In this case, the operator T was a FIO with phase 
function of quadratic type, cf. p], having natural action on the Gabor atoms g m ,n- 

The treatment of hypebolic problems in Gabor frames remained outside the results 
of |14j . In this regards, the example of the transport equation with non-constant 
coefficients, standard test for numerical schemes, is somewhat discouraging. In fact, 
for a fixed time t > 0, the operator T solving the Cauchy problem turns out to be 
a change of variables, and Gabor frames are not suited to follow the corresponding 
non-linear propagation of singularities. 

Because of this, our attack to hyperbolic equation via Gabor frames will be limited 
to operators with constant coefficients. Nevertheless, with respect to [3 [6], our results 
will be stronger in the following aspects: 

• we shall treat hyperbolic equations, not necessarily strictly hyperbolic, of any 
order and in an arbitrary number of space variables; 

• our class of equation will include also parabolic equations; 

• the decay of the Gabor matrix will be of super-exponential type. 

As a simple example, let us consider the Cauchy problem for the wave equation: 

(4) dfu - A x u = 0, (t, x) g R x R d , 

it(0, x) = uo(x), dtu(0,x) = u\(x). 

We may express the solution in the form 

(5) u(t, x) = T t u\(x) + d t T t uo(x), 
where Tt is the Fourier multiplier 

(6) T t f(x) = J e 2 ^a t (0f(0di, 
with symbol 

(T) *(«)- 8h », {6*. 

Taking the Gaussian as window function g and discretizing consequently T = Tt, we 
get a Gabor matrix satisfying 

(o\ IT / / I < f7 p -e(|m-m'| 2 + |n-n'| 2 ) 

V°7 \- L m ,n ' ,m,n\ — 
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We address to the sequel of the paper for the dependence of C and e on t > 0. With 
respect to the approach by FIOs, the ^-dependent phase functions are here absorbed 
into the symbol ([7]). Because of its oscillations, 0t(£) then belongs the the non-standard 
Hormander's class Sn , defined in general by the estimates 



0) 



\d a a{z)\ < C a 



(x,0 G 



i,2./ 



As well known, PSDOs with symbol in Sq have unpredictable action on singular 
supports. In spite of this, the estimates are very strong. 

Let us describe in short how ^ can be obtained in the general context of symbols 
depending on z = The proof will be largely based on [221 EU [32], where the 

authors show that for every PSDO with symbol a G S® , taking as window any function 
g in the Schwartz space we have 



(10) 



m' ,rr ,m,n 



< C N (1 + \m 



'|2 i I 

m + In 



n2\-N 



for N as large as we want. Here, inspired by the properties of the symbol in ([7]) we 
introduce an analytic version of the class Sqq, by assuming for s > 



(11) 



\d a a{z)\ < C^ +l {a\) 



z G 



D 2J 



for a constant C > depending on a. We say that the symbol a is Gevrey when s > 1, 
analytic when s = 1, and ultra-analytic if s < 1. As we will show, such estimates 
are satisfied by the symbols of the fundamental solutions for a large class of evolution 
operators with constant coefficients, with < s < 1. In particular, we have s = 
for the wave equation Q, and s = 1 — l/(2k), k G N \ {0}, for the generalized heat 
equation: 



(12) 



d t u + (-A x ) k u = 0, u(0,x) = u o . 



Let T be a PSDO with symbol a satisfying (11) for some s > 0, and define 



min{2, 1/s}. Fix as a window the Gaussian, or more generally any function g belonging 
to the Gelfand-Shilov classes S^.(R. d ). Then we shall prove that 



(13) 



-m' ,n' ,m,n 



< Ce 



t{\m—m'\+\n—n'\) T 



Accordingly, every column or row of the Gabor matrix, rearranged in decreasing order, 
displays a similar decay, i.e. we obtain an exponential-type sparsity. 



In the case of the Gevrey-analytic symbols, i.e. s > 1 in (11) and r < 1 in (13), the 



estimates follow from the results in [23]. In view of the applications, our attention will 
be mainly addressed to the ultra-analytic case, i.e. < s < 1 in (11), r > 1 in (13); in 



particular from ( 13 ) with r = 2 we have (1ST) for the wave equation. 



The plan of the paper is the following. In Section 2 we provide the necessary prelim- 
inaries, concerning Gelfand-Shilov spaces, time-frequency representations and modula- 



tion spaces. Section 3 is devoted to the study of the symbols in ( 11 ). The corresponding 



PSDOs are considered in Section 4, where we prove ( 13 ) and give boundedness results 
in modulation and Gelfand-Shilov spaces. Applications to evolution equations are 
given in Section 5. 
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2. Preliminaries 

2.1. Notations. The Schwartz class is denoted by the space of tempered distri- 

butions by S'(R d ). We use the brackets (/, g) to denote the extension to S'(R d ) x S(R d ) 
of the inner product (f,g) = J f(t)g(t)dt on L 2 (R d ). 

Euclidean norm of x G R d is given by \x\ = (x\ + • • • + x 2 ^ 1 ^ 2 , and (x) = (1 + 
Ixj 2 ) 1 / 2 . We write xy = x ■ y for the scalar product on R d , for x, y G R d . 

For multiindices a, /3 G N d , we have |a| = cci + • • • + a^, a! = ai! • • -a^!, x a = 

^-x^.xe M d , and, for /3 < a, i.e. & < a„ j € {1, 2, ... , d}, (») = (£) (^) . 

The letter C denotes a positive constant, not necessarily the same at every appearance. 

The operator of partial differentiation d is given by 

f) a = f) a = f) ai ■ ■ ■ f) ad 

for all multiindices a £ N d and all x = (x±, . . . , x^) G If / and g are smooth enough, 
then the Leibniz formula holds 



d a (fg)(x) = Y, (fyd"-ef(x)dPg(x). 



/3<a 

Observe 
and 

(15) |a|! < dl Q la!. 

The Fourier transform is normalized to be /(£) = Ff{£) = f f(t)e~ 2mt ^dt. 
Translation and modulation operators, T and M are defined by 

T x f(-) = f(--x) and M x f(-)=e 27rix -f(-), x G M d . 

For z = (x,£) we shall also write 

ir(z)f = M(T x f. 

The following relations hold 

M y T x = e 2nixy T x M y , (T x fJ=M. x f, (M x ff=T x f, x, y G R d , f, g G L 2 (R d ). 

For < p < oo and a weight m, the space -d(A) is the space of sequences a 
{«a}asA on a lattice A, such that 

/ \ Vp 



E KIMA) P 

vAgA y 



< OO 



(with obvious changes when p = oo). 

Throughout the paper, we shall use the notation A < B to express the inequality 
A < cB for a suitable constant c > 0, and A x B for the equivalence c~ 1 B < A < cB. 
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2.2. Gelfand-Shilov Spaces. The Schwartz class S(M. ) does not give enough infor- 
mation about how fast a function / G S(M. d ) and its derivatives decay at infinity. This 
is the main motivation to use subspaces of the Schwartz class, so-called Gelfand-Shilov 
type spaces, introduced in |20j and now available also in a textbook |30j . Let us recall 
their definition and main properties. 

Definition 2.1. Let there be given s,r > 0. A function f G is in the Gelfand- 

Shilov type space S^(M. d ) if there exist constants A,B>0 such that 

(16) \x a dPf(x)\ <A\ a \B^(alY(p\y, a,/3GN d . 

The space S^(R d ) is nontrivial if and only if r + s > 1, or r + s = 1 and r, s > |20j . 
So the smallest nontrivial space with r = s is provided by S^^O^)- Every function of 

the type P(x)e~ a ^ 2 , with a > and P(x) polynomial on R d , is in the class S^,^(M. d ). 
We observe the trivial inclusions S^(R d ) C S^(R d ) for s\ < S2 and r± < r-i- Moreover, 
if / G S^(R d ), also x s d^f belongs to the same space for every fixed 5,j. 

The action of the Fourier transform on S^.(R d ) interchanges the indices s and r, as 
explained in the following theorem. 

Theorem 2.2. For f G S(R d ) we have f G S°(R d ) if and only iff G S r s (R d ). 

Therefore for s = r the spaces Sg(R d ) are invariant under the action of the Fourier 
transform. 

Theorem 2.3. Assume s > 0,r > 0, s+r > 1. For f G <S(M ), the following conditions 
are equivalent: 

a) / G S s r (R d ) . 

b) There exist constants A,B > 0, such that 

\\x a f\\ L °o<AW(a\y and \\^f\\ L oo < B^(/3\) s , a,/3GN d . 

c) There exist constants A, B > 0, such that 

\\x a f\\ L ~<AM(a\y and ||^/|| L oc < BW(J3\)', a,/?GN d . 

d) There exist constants h, k > 0, suc/i f/iai 

||/e /l N 1/r || £;00 < oo anrf ||/ e fc l€l 1/s || L oo < oo. 

A s uita ble window class for weighted modulation spaces (see the subsequent Defi- 
nition 2.6) is the Gelfand-Shilov type space S|(M rf ), consisting of functions / G 5(M d ) 



such that for every constant A > and B > 

(17) \x a d p f{x)\ < AMBWalpl, a,/3e N d . 

We have S s JR d ) C X^IR ') C 5j(M d ) for every s < 1. Observe that the characterization 
of Theorem 2.3 can be adapted to S}(M d ) by replacing the words "there exist" by "for 
every" and taking r = s = 1. 

Let us underline the following property, which exhibits two equivalent ways of ex- 
pressing the decay of a continuous function / on R d . This follows immediately from 
[301 Proposition 6.1.5], where the property was formulated for / G S(W i ). For sake of 
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clarity, we shall detail the proof showing the mutual dependence between the constants 
e and C below. 

Proposition 2.4 ([30, Proposition 6.1.5]). Consider r > and let h be a continuous 
function on R d . Then the following conditions are equivalent: 

(i) There exists a constant e > such that 

(18) \h(x)\ < e- e ^\ x£R d , 

(ii) There exists a constant C > such that 

(19) \x a h{x)\ < Cl Q l(a!) r , x£R d ,a£N d . 



Proof. We re- write ( 18 ) in the form 
(20) 



\h(x)\r < e-rM*, xGK°\ 



In turn, (20) can be re- written as 

oo 

(21) 

Hence the sequence of the terms of the series is uniformly bounded, as well as the 
sequence of the r-th powers: 



sup (-^) (n!) 1 \x\ r \h(x)\ r < oo. 



-(n!)- r |xn/i(x)|, nGN, 



and we obtain 



Hence, writing \a\ = n and applying (|15j) : 

' j\ Hal 
ra x 1 1 



\x\ n \h(x)\<—(niy, xe 



n = 0, 1,.... 



\x a h(x)\ 



(a\) r = C\ a \{a\) r , x£R d , a G N d 



where C 



Therefore (19) is proved. 



In the opposite direction, let (19) be satisfied. Using the following inequalities 



M"<Esi*"i. E 



a. 



d n , and a! < \a\l, 



\a\=n \a\=n 

and the assumption |l9| ) we obtain 

bri/i(x)|< V ^|x Q /i(rr)|< V ^C H (a!) r < C n {n\) r d n = (dC) n (n\) r , 

\a\=n \a\=n 

for every x G n G N. Therefore the sequence 

(dC)-"(n!)- r |2;| n |/i(x)|, n G N 
is uniformly bounded for x G as well as the sequence 

(dC)~r(n!)- 1 |x|?|/ l (x)|?, nGN. 
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If we choose e = q{dC)~r ; for a fixed q G (0, 1), we conclude 

j oo oo 

e eN? |/i(x)|r = ^(dCr^n!)- 1 !*!^)!' < 



n=0 n=0 



This is (20), hence (18) and the proof is complete. □ 



It follows from this proof the precise relation between the constants e and C. Indeed, 



assuming (18), then (19) is satisfied with C = I — I . Viceversa, (19) implies (18) for 



any e < r(dC) r . The bound is sharp for d = 1. Also, it follows from the proof that 
the constant implicit in the notation < in (|18|) depends only on the corresponding one 



in ( 19 ) and viceversa. 

The strong dual spaces of and X];(IR rf ) are called spaces of tempered ultra- 

distributions and denoted by (S^)'(M. ) and respectively. Notice that they 

contain the space of tempered distribution S'(M. ). 

The spaces S^(R d ) are nuclear spaces [29]. This property yields a kernel theorem 
for Gelfand-Shilov spaces, cf. [37J. 

Theorem 2.5. There exists an isomorphism between the space of linear continuous 
maps T from S%(M. ) to (S^)'(R d ) and (S^:)' (R 2d ) , which associates to every T a kernel 
K T G {Sf)'{R 2d ) such that 

(Tu,v) = (K T ,v®u), Vu,v £ S s r (R d ). 

Kt is called the kernel ofT. 

2.3. Time- frequency representations. We recall the basic concepts of time-frequency 
analysis and refer the reader to [21] for the full details. Consider a distribution 
/ G e>'(IR d ) and a Schwartz function g G 5(M a! ) \ {0} (the so-called window). The 
short-time Fourier transform (STFT) of / with respect to g is V g f(z) = (f,ir(z)g), 
z = (x, £) G K M . The short-time Fourier transform is well-defined whenever the 
bracket (•,•) makes sense for dual pairs of function or (ultra-) distribution spaces, in 
particular for / G S'(R d ) and g G S(R d ), f,g€ L 2 (R d ), f G (S})'(M d ) and g G S}(M d ) 
or / G (5y)'(]R ) and g G S^(M. d ). Let us recall the covariance formula for the short- 
time Fourier transform that will be used in the sequel 

(22) V g (M v T y f)(x, = e~ 2m ^-^(V g f)(x - y,S - rj), x, y,£,T) G R d . 

Another time-frequency representation we shall use is the (cross-)Wigner distribu- 
tion of f,g G L 2 (R d ), defined as 



(23) W(f,g)(x, t)=f d f(*+l)g(x- \)e~ 2 ^ dt. 

If we set g{t) = g(—t), then the relation between cross- Wigner distribution and short- 
time Fourier transform is provided by 

(24) W(f, g)(x, = 2 d e^V s f(2x, 2$. 

For the discrete description of function spaces and operators we use Gabor frames. 
Let A = AL 2d with A G GL(2d, R) be a lattice of the time-frequency plane. The set 
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of time-frequency shifts Q(g, A) = {-7r(A)g : A G A} for a non-zero g G L 2 (K d ) is called 
a Gabor system. The set Q(g, A) is a Gabor frame, if there exist constants A, B > 
such that 

(25) < £ |(/,vr(A) 5 )| 2 < B||/||2, V/ G L 2 (R d ). 

AeA 

If (25) is satisfied, then there exists a dual window 7 G L 2 (M a! ), such that Q(j, A) is a 
frame, and every / G L 2 (M ci ) possesses the frame expansions 

/ = ^(/,vr(A)<7)7r(A)7 = £</, vr(A) 7 )vr(A) 5 
AeA AeA 

with unconditional convergence in L 2 (IR d ). In conclusion, we list some results about 
time-frequency analysis of Gelfand-Shilov functions, cf. [8l \15 \ [251135] : 

(26) f,ge S s s (R d ), s > 1/2 V g f G Sf(R 2d ), 

(27) f,ge S s s (R d ), s > 1/2 G S s s (R 2rf ), 

(28) /, 5 G Si (M d ) W (/, 5 ) G Sj (R 2d ) . 
If p G S s s {R d ), s > 1/2, then 

(29) / G 5f(M d ) |V fl (/)(*)| < e" e|2|1/s for some e > 0. 

2.4. Modulation Spaces. Weighted modulation spaces measure the decay of the 
STFT on the time-frequency (phase space) plane and were introduced by Feichtinger 
in the 80 's |19| . for weight of sub-exponential growth at infinity. The study of weights 
of exponential growth at infinity was developed in |15} 135] , so that we can now recall 
their definition for weights with at most exponential decay at infinity. 

Weight Functions. In the sequel v will always be a continuous, positive, even, 
submultiplicative function (submultiplicative weight), i.e., u(0) = 1, v(z) = v(—z), 
and v{z\ + z%) < v(zi)v(z2), for all z,z\,Z2 G M. 2d . Submultiplicativity implies that 
v(z) is dominated by an exponential function, i.e. 

(30) 3C, k > such that 1 < v{z) < Ce fc|z| , z G M. 2d . 
For example, every weight of the form 

(31) v(z) = e s ^\l + \z\) a log r (e + \z\) 

for parameters a,r, s > 0, < 6 < 1 satisfies the above conditions. 

We denote by A4 v (K d ) the space of u-moderate weights on WL d ; these are measurable 
positive functions m satisfying m(z + C) < Cv{z)m{Q for every z,( G M. d . 

Definition 2.6. Given g G S}(M d ), a weight function m G M. v (J$? d ), and 1 < p, q < 
00, the modulation space Mm 9 (^ d ) consists of all tempered ultra- distributions f G 
(S^)'(]R d ) such that V g f G L%?(R 2d ) (weighted mixed -norm spaces). The norm on 
M^ q (R d ) is 



(32) 11/11^ = 11^/11^=1/ (/ \V g f(x,0\ p m(x,0 P dx) <%) 
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(obvious changes if p = oo or q = oo). 

We observe that for /,g £ Sj;(]R d ) the above integral is converg ent and thus £} (R d ) C 
Mm Q (R d ), 1 < p, q < oo, cf. [IB] , with dense inclusion when p,q < oo, cf. [8]. When 
p = q, we simply write M^(M d ) instead of M%l p (R d ). The spaces M% q (R d ) are Ba- 
nach spaces and every nonzero g G M^K^) yields an equivalent norm in ( |32[ ) and so 
M%, g (R d ) is independent on the choice of g G M^(R d ). 

We recall the inversion formula for the STFT (see ([HJ Proposition 11.3.2] and [151 
Proposition 2.6] for exponential weights): assume g G M^(R d ) \ {0}, / G Mm Q (R d ), 
then 

(33) /=_L/" V g f(x,t)MtT x gdxd£, 

\\9h J^ 2d 

and the equality holds in M^ q (R d ) (observe that M^(R d ) C M 2 (R d ) = L 2 (R d ), so 
\\g\h < oo). 



3. Gevrey-analytic and ULTRA-ANALYTIC SYMBOLS 

In this section we characterize the smoothness and the growth of a function / on 
R d in terms of the decay of its STFT V g f, for a suitable window g. In the proofs we 
shall detail the relations among the constants Cf, C g , Cf >g and e which come into play, 
since it could be useful for constructing efficient algorithms related to these objects. 

Theorem 3.1. Consider s > 0, m G M v (R d ), g G M^^W 1 ) \ {0} such that there 
exists C g > 0, 

(34) II^IIlJ(r*)< CW(a!)», a G N d . 

For f G C°°(R d ) the following conditions are equivalent: 

(i) There exists a constant Cf > such that 

(35) \d a f(x)\ < m(x)c} a| (a!) s , x G R d , a G N d . 

(ii) There exists a constant Cf g > such that 

(36) \eV g f(x,0\ <m(x)C ] ^ g (aV) s , (x, £) G R 2d , a G N d . 
(Hi) There exists a constant e > such that 

(37) \V g f(x,0\<m(x)e-^ h , (x, £) G M M , a G N d . 



if i/ie equivalent conditions (|35l), (|36|), (37) are satisfied, we say that f is a Gevrey 



symbol when s > 1, an analytic symbol if s = 1 and an ultra- analytic symbol when 
s < 1. 

Proof, (i) =>• (ii). Since f,T x g G C°°(M rf ), for every x G M d , we can use Leibniz formula 
and write 



d a (fT x g)(t) = ^ ( a \d a -Pf{t)T x (d?g)(t), teR d ,aen d 

/3<a 
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Let us estimate the L 1 (M' i )-norm of the products d a ~^ fT x (d@g). Using the positivity 
and u-moderateness of the weight m, Holder inequality and the assumptions (35) and 
©, 

\\ff*- /, fT m (a' , g)\\ IA < \\d^f\\ LT/ JmT x (d^g)\\ Ll 

< C f- p \(a - P)\) S \H- - x)dPg{. - x)\\ L im(x) 

< C^(a\) s m(x) 

where C := max{Cf,C g } and we used (a — 0)\j3\ < a!. These estimates tell us, 
in particular, that V g f(x,£) = F{fT x g){£) is well-defined. We can exchange partial 
derivatives and Fourier transform as follows 



Using 



\T(d a -PfT x (dPg))(0\ < \\Hd a - p fT x {dPg))\\ L ~ < \\d a ^ fT x (d^g)\\^, 
the majorizations above and ( |14[ ), 

k-v./(,,oi < (fjc^'-r = ^cH(d)-. 



max{Cy, C 9 } 



7T 



This proves (36), with constant C/ j9 : = 

(ii) (i). We use the inversion formula (33), observing that (36) implies / € 
M^ 1 (M c( ), so that the equality in (33) holds a.e. and we can assume that it holds 
everywhere since / is smooth. Let us consider the partial derivatives of / and ex- 
change them with the integrals, the estimates below will provide the justification of 
this operation. So, formally, we can write 



d a f(t) 



1 



\g\\i 



V g f(x,Od a (M i T x g)(t)dxdC, t€ 



Using Leibniz formula d a {M{T x g) = Y,p< a © (2iri^ M^id 01 '! 3 g) we estimate 

\d a f(t)\ < A E («V 2vr ) l/3 ' / l^/M^I ■\T x d a ^g(t)\dxdi ) tGM 
\\9h \PJ J^ 2d 



We set 



Ta,p(t) := / \V g f{x,0^\ ■ \T x d a -Pg(t)\ dxdi 

jR 2d 



and prove that, for every fixed t G I a 8(t) are absolutely convergent integrals. 
Using 1 + < CiEm <d+i l£ 7 l> where Cd depends only on the dimension d, and 
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assumption (36) 

la A*) 



1 _|_ 

m \V g f(x,Of\ g \Txd^- p g{t)\ dxdi 



|t|<<Z+1 



| 7 |<d+i 



Now, 09 + 7)! < 2l^l+lTl( ( 0!)( 7 !) < 2 d+1 {d + 1)12^01 and cf+ jl = Cffcf^ so, forget- 
ting about the constants depending only on d, f and g, 

IaAV £ Cf^iPiy (J m{x) \T x d a -Pg{t)\ dx 

< Cf^((3l) s m(t) U v(x - t)\d a -?g(t -x)\dx 
= C^ g 2 s P{p\) s m{t)\\d a -P 



9Ul> 



The estimate above and the assumption (34) on g allow the following majorization: 



\d a fit)\ 



< 



< 



< 



m(t) 

Nit 

m(t) 



/3<a 



E U (2t)'^' (2 s C /i9 )I^I (/3!)«Cj-fl ((a -0)!)' 



(a!) s max{C /i9 ,C 9 }l a l2 s H(27r) |a| 2H 



m(i)C/(a!) s , 



where C/ := 2 s+i 7rmax{C/ i9 , C g }, and we used (/3)!(a — /3)! < a! and (14). 
(u) 44> (m). The equivalence is an immediate consequence of Proposition 2.4 and 
the subsequent remarks on the constants, where, for every fixed x G M. d , we choose 
h(£) := V g f(x, £)/m(x) and r = s. □ 



A natural question is whether we may find window functions satisfying (34). To 
this end, we recall the following characterization of Gelfand-Shilov spaces. 

Proposition 3.2. Let g G S(R d ). We have g G S?(R d ), with s,r > 0, r + s > 1, if 
and only if there exist constants A > 0, e > such that 

\d a g{x)\ < A^{al) s e- e ^\ x G R d , a G N d . 

We have g G if and only if, for every A > 0, e > 0, 

\d a g(x)\ < A^ale' 6 ^, x G M. d , a G N d . 

Proof. The first part of the statement is in [30} Proposition 6.1.7]. For the second part, 
the assumption g G Sj;(R ) means that for every A > 0,B > 0, 

\x p d a g{x)\ < A^B^aipi, x£R d , a, /3 G N d . 
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Therefore the function h 



d a g 
A\ a \al 



satisfies (19) in Proposition 



2.4 



for every C > 0, 

and the estimate (18) is then satisfied for every e > (see the observations after the 



proof of Proposition |2.4| for the uniformity of the constants). This gives the claim. □ 

""^ with s > 0, < r < 1, s + r > 1, satisfies (34) for every 



Hence every g G Sf,(I 
submultiplicative weight v (see (30)). The same holds true if g G S^(M d ) and s > 1. 



4. Almost Diagonalization for Pseudodifferential operators 

In this section we extend some results of [21] to the case of pseudodifferential oper- 
ators having (ultra-) analytic symbols. Our result covers also the Gevrey-analytic case 
which was already discussed in |24j . 

The Weyl form a w of a pseudodifferential operator (so-called Weyl operator or Weyl 
transform) with symbol cr(x, £) on R 2d is formally defined by 



a 



y 



ri)e 2 *« x -y^f(y)dyd V . 



Using the Wigner distribution (23), we can recast the definition as follows 
(39) (a w f,g) = (a,W(f,g)). 

Since for f,g€ S{R d ) we have W(f,g) G S{R 2d ), then ((39]) defines for a G S'(R 2d ) a 



continuous map o~ w : 5(M a! ) — > S'(R d ). The Schwartz kernel K of a w is given by 



(40) 



K(x,y) 



T x \ x _ y a 



,£) ^S\R 2d ) 



On the other hand, in view of the kernel theorem in S — S' , every linear continuous 
map from 5(M d ) to 5'(IR M ) can be represented by means of a kernel K G 5'(M 2d ), 



hence by means of (38), (39) with symbol 
(41) 



a(x, £) 



F y -+t-K[x + -,x 



Same arguments are valid in Gelfand-Shilov classes with s = r > 1/2 (see |29^ I37j). 
Namely, considering a G {S s s )'{R 2d ), f,g£ S s s {R d ) in (g, ((39}, we have a w : S s s (R d ) 



(Sg)'(R ) continously, in view of (27). The kernel K of it " 1 , g iven by (40), belongs to 
(Sg)'(R 2d ). In the opposite direction, in view of Theorem 2.5, every linear continuous 
map from 5f (R d ) to (S s s )'(R d ) can be represented in the form @, ([39]), with a G 
(Sg)'(R 2d ) given by (41). The same holds for the couple of spaces Ti\, (S^)'. 



The crucial relation between the action of the Weyl operator a w on time-frequency 
shifts and the short-time Fourier transform of its symbol, contained in [221 Lemma 3.1] 
can now be extended to Gelfand-Shilov spaces and their dual spaces as follows. 

Lemma 4.1. Consider s > 1/2, g G S s s (R d ), $ = W(g,g). Then, for a G {S s s )'{R 2d ), 



(42) 
and 
(43) 



\(a w 7T(z)g,7r(w)g)\ 



V$CT 



| V$tr(n, v ) 



(a w ir u 



z + w 



g^ 



\V$(t(u,v)\ 



u H — j~ 
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where j(zi, Z2) = (z2, —z\). Moreover, the same results hold true if we replace the space 
S s s (R d ) with the space S^(]R d ). 

Proof. Since $ = W(g,g) G S s s (R 2d ) for g G S s s (R d ) the duality (a, tt(u, w)$) (S| y xS j 
is well-defined so that the short-time Fourier transform V$o~(u,v) makes sense. The 
same pattern applies to the case g G X^(IR d ). The rest of the proof is analogous to [221 
Lemma 3.1]. □ 

We now exhibit a characterization of smooth symbols a G C OD (R 2d ) for o~ w in terms 
of the continuous Gabor matrix. The symbol a belongs to the classes defined in Section 
3, the dimension being now 2d. 

Theorem 4.2. Let s > 1/2, and m G M. v {R 2d ). If 1/2 < s < 1 consider a window 
function g G Sg(R d ); otherwise, if s > 1, assume either g G or g G Sg(R d ) 

and the following growth condition on the weight v: 



,{z)<e^'\ zG 



i,2tl 



for every e > 0. Then the following properties are equivalent for a G C QC (R 2d ): 
(i) The symbol a satisfies 

(44) \d a a(z)\<m(z)C M (a\) s , Vz 6 R 2d , Va G N 2d . 

(ii) There exists e > such that 



(45) 



\(a w Tr(z)g,ir(w)g)\ < m 



w + z 



3 — e\w — z\ 1 



V z, w G 



1)2(1 



Proof, (i) (ii). Proposition 3.2 applied to the window in Lemma 4.1 and the 



assumptions on v (recall also (30)) imply that <3? satisfies the assumptions of Theorem 
3.1 Hence, using the equivalence p5|) 44> (37), the assumption (44) is equivalent to 



the following decay estimate of the corresponding short-time Fourier transform 



\V*a(u,v)\ < m{u)e- e ^\ u,v G 



i,2d 



for a suitable e > 0, which combined with (42) yields 

,z + w 



\ V ^ a {-^—^{w - z))\ < m 



w + z 



-e\j(w-z)\s 



m 



w + z 



-e\w — z\ 1 



that is (ii) 



(ii) (i). Relation (43) and the decay assumption (45) give 



I V$cr(n, v)\ 



(o- w ir(u-^j \v)\g,7r(u + ^j \v)\ g) 



< 



m(u)e 



m(u)e 



-e|t)|1 



and using the equivalence (35) 44> (37) we obtain the claim. 



□ 



Of course, from (45) we deduce the discrete Gabor matrix decay in (48) below. 



The viceversa requires more effort. Indeed, we appeal to a recent result obtained 
by Grochenig and Lyubarskii in |23j . They find sufficient conditions on the lattice 
A = AT?, A G GL(2,M), such that g = Y^k=o c kHk, with Hermite function, forms a 
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so-called Gabor (super)frame Q(g, A). Besides they prove the existence of dual windows 
7 that belong to the space Sy 2 (M.) (cf. \23\ Lemma 4.4]). This theory transfers to the 

1 /2 rf 

d-dimensional case by taking a tensor product g = gi <S> • • • <S> <?d G Sy 2 (Mr) of windows 
as above, which defines a Gabor frame on the lattice Ai X • • • x and possesses a 
dual window 7 = 71 (g) • • • (g) 7^ in the same space G sH* Let us simply call Gabor 
super-frame Q(g,A) for L 2 (M d ) a Gabor frame with the above properties. The Gabor 



super- frames are the key for the discretization of the kernel in (45). First, we need the 
preliminary result below, which reflects the algebra property of Gelfand-Shilov spaces. 
For s > 1/2, e > 0, we define the following weight functions: 



(46) w 8>e (z) :=e-^ s , zG 



Lemma 4.3. Let A C M? d be a lattice ofM? d . Then the sampling {w s ^(X)}x e /^ of (46) 
( defined on M. 2d ) satisfies 

, . . , , , v-^ /^ \ 1 \ ^ [ w s e(A), for s > 1 

(47) K e *^, e )(A):=^^, e (A-z/K>)<| ^ £2 _ i/a(A)) ^ i< fi<1 

Proof. For s > 1 the result immediately follows by the submultiplicativity of w s<e : 
(w Sje * u> s , e )(A) = ^ i« s ,e(A - v)w s ^{v) < io S)£ (A) ^ w Si£ (-^)u; Sie (z;) < w a , e (A). 

For 1/2 < s < 1, we use the arguments of |21l Lemma 11.1.1(c)]. For A G A, we divide 
the lattice A into the subsets N\ = {v £ k : |A — v\ < |A|/2} and iVjj = {z/ G A : 
|A - u\ > |A|/2}. For v G AT A , \v\ > |A|/2 and ji/j 1 / 5 > (IAI/2) 1 / 5 , so 

K, e *^, e )(A)<e-( e2 " i )l A l" ( J2 e-^° + e- eH M < e~^~'^' . 

This concludes the proof. □ 

Theorem 4.4. Let Q(g, A) a Gabor super-frame for L 2 (M. d ). Consider m G M. v (Jst 2d ), 
s > 1/2, and a symbol a G C°°(R 2rf ). T/ien t/ie following properties are equivalent: 



(i) There exists e > such that the estimate (45) holds. 

(ii) There exists e > such that 



(48) |(er w 7r(/i)#,7r(A)#)| < m e~^~^ , VA,/i G A. 

Proof. It remains to show that (ii) =^> (i). The pattern of |22| Theorem 3.2] can 
be adapted to this proof by using a Gabor super-frame Q(g, A), with a dual window 

7 GS 1 1 / / 2 V). 

Let Q be a symmetric relatively compact fundamental domain of the lattice A C M. 2d . 
Given w,z G M. 2d , we can write them uniquely as w = A + u, z = fi + u', for A, ji G A 



1 



GABOR REPRESENTATIONS OF EVOLUTION OPERATORS 



15 



and u, v! G Q. Using the Gabor reproducing formula for the time-frequency shift 
n(u)g G Sy^ 0^ d ) we can write 

i/EA 

Inserting the prior expansions in the assumption ( 48 ) , 

\{a w ^ + u')g^{\ + u)g)\ 

< K^vr(/i + ^^(A + ^)ll(vr(n05,vr(^7)IK^)5^H7)l 

u,u'eA 

(49) < m( X + fl+ 2 U + U ' ) e-^-^'W^ - u>)\\V 7 g(v - u)\. 



v,v'&A 



Since the window functions g, 7 are both in Sy^(M. d ), the STFT V^g is in S^%(M?"*) in 

' ■ '-' ■■■hi 



l>2d\ 



view of (26). Thus there exists h > such that |V^g(^)| < e _/l ' 2 ' , for every z G 
In particular, being Q relatively compact and u G Q, 

l^ 7 9(^ "-")!< e^ 1 ^ 2 < sup e-^'l 1 -"! 2 < e - h M 
The assumption m G M v (R 2d ) yields 



u6Q 



X + fi + u + u'\ / X + u\ (v\ ( v' 
< m 1.1 



w — u 



and, for every < h < h, 



Inserting these estimates in (49), 



-%| a < e -ftM a Vi/gA. 



|(a to 7r( M + «')<?, 7r(A + u) 5 }| 



< 



m 



2 



< . i' A • /' 



(50) 



m 



< m 



' u,u'eA 

\ Y e -b(\X+u-^,--u'\l+\iy\ X s+\u'\l) 



e -m\ ( V _ \ p -h\»'\ 2 



with b = e for s > 1/2 whereas b = min{e, h} for s = 1/2 (there may be a loss of 
decay) . We observe that the raw ( 50 ) can be rewritten as 



m 



( w s,b * w s,b * w Stb )(X - fi), 
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-£>|A| 1/s _ Njow we apply Lemma 



4.3 



twice and we obtain 



(50) <m 



-e\\-n\s 



with e = b when s > 1 whereas e = 62~ 2 / s when 1/2 < s < 1. Again, observe the loss 
of decay for 1/2 < s < 1. 

If w, z G M M and w = \+u, z = fi+u', A, \i G A, u, it' G Q, then A — \i = w — z+u' — u 
and u' — u G Q — Q, which is a relatively compact set, thus 



r - '• - sup e -E|t«— 1 < e -e\w-z\l 



(52) e -£]A-^[. < 

«eQ-Q 

Finally, the f -moderateness of the weight m G 7Wt)(IR 2c( ), together with the fact that v 
is continuous and the set Q + Q is relatively compact, let us write 



(53) 



771 



A + V 



771 



< 



777 



W + Z u + v! 
w + z 



< 



u 



w + z 



u + u' 



SU P v v ■ > 

ueQ+Q KZ 



< 



m 



w 



Combining the estimates (52) and (53) with (51) we obtain (45), with the parameter 
e = i which appears in (51 ). □ 

4.1. Sparsity of the Gabor matrix. The operators a w which satisfy Theorem |4.2[ 
say with m = v = 1, enjoy a fundamental sparsity property. Indeed, let Q(g,A) be a 
Gabor frame for L 2 (R d ), with g £ S s s (R d ), s > 1/2. Then, as we saw, 

(54) | (cr™ 7r(/j)g, 7r(A)g)| < Ce" 6 !^', V A, /i G A, 

with suitable constants C > 0, e > 0. This gives at once an exponential-type sparsity, 
in the sense precised by the following proposition (cf. [26] for the more standard 
notion of super-polynomial sparsity). 



Proposition 4.5. Let the Gabor matrix (a w ir({i)g,7r(\)g) satisfy (54). Then it is 
sparse in the following sense. Let a be any column or raw of the matrix, and let \a\ n 
be the n-largest entry of the sequence a. Then, \a\ n satisfies 

\a\ n <Ce- en ^, 7i GN 
for some constants C > 0, e > 0. 
Proof. By a discrete analog of Proposition |2.4| it suffices to prove that 

n a \a\ n < C a+1 (al) 2ds , a G N. 

On the other hand we have 

i 

nv ' \o\n ^ ll a ll^) 

for every < p < oo. Hence by (|54[) and setting p = 1/a we obtain 



n a \a\ n < ( ^ e 



-ep\\-n\ 



AeA 



AeA 



-ep|A|a\ v 
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Let Q be a fundamental domain of the lattice A. Then if x G A + Q, A G A, we have 
\x\ < |A| + C , therefore (x] 1 / 3 < CidA] 1 / 3 + 1). Hence 

^2 & -e P \\\s <Q 2 J e -ep\x\s dx= j da I e - epps p2 d-l df) 



AeA 

-+oo 



C 3 s [ + ™ e -t t 2ds-i dt _ C 3 sT(2ds) _ C 4 



(ep) 2ds J (ep) 2ds p 2ds 

Finally, by Stirling's formula, 

n a \a\ n < < C« +1 (a!) 2ds . 

p p 



□ 



4.2. Boundedness of pseudodifferential operators on modulation and Gelfand- 
Shilov spaces. Although not necessary for our numerical applications, we list some 
continuity results on modulation and Gelfand-Shilov spaces which follow easily from 
our basic Theorem 021 



Proposition 4.6. Let s > 1/2 and consider a symbol a G C°°(M. 2d ) satisfying the 
estimates 

(55) \d a a(z)\<C^(al) s , 

for some C > 0. Let m G A4 v (M? d ) and, if s > 1, assume the weight v satisfies 

v(z)<e^ 1/a , z G M. 2d 

for every e > 0. Then the Weyl operator a w extends to a bounded operator on 
M™(R d ). 

Proof. Let g G £{(1^) with ||g||£2 = 1. From the inversion formula (33), 



V 9 (a w f)(u) = / (a w n(z)g, ir(u)g) V g f{z) dz. 
The desired result thus follows if we can prove that the map M{a) defined by 
M{a)G{u)= [ (a w Tr(z)g,7r(u)g)G(z)dz 



4.2 



is continuous from Lm 9 (K 2a! ) into L%f(R 2d ). The characterization of Theorem 
assures the existence of an e > 0, such that 

\V g (a w f)\ = \M(a)V g (f)\ < (w., e * \V g f\)(u), u G R 2d 

where the weight function w Sj6 on M? d is defined in ( |46[ ). The desired conclusion then 
follows from the relation L V m * L^ <-} L™ , for w Sj6 G Lr v by the growth assumption on 
v (if s < 1, (30) suffices). 

□ 

Proposition 4.7. Let s > 1/2, and consider a symbol a G C°°(M. 2d ) that satisfies 
(55). Then the Weyl operator a w is bounded on 5|(IR d ). 
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Proof. Fix a window function g G (R d ). For / G S s s (R d ) we have V g f G S s s (R 2d ) 



in view of (26). Hence there exists a constant h > such that < e ^l 2 ' " , for 



eve ry z G M . Thus, taking e = min{e, h} and using a continuous version of Lemma 
4.3 for every u G R 2d , 

\V g (a w f)(u)\ < (w s , € *w Sih )(u) < (ws^ws^iu) - t "'" A " 1 ' 0r * - ' 



iu s gj-VsC^)) f° r | < s < 1. 



Using (29) we obtain the claim. □ 



Boundedness results for pseudodifferential operators on Gelfand-Shilov spaces are 
contained in |36j. 

5. Applications to evolution equations 

In this section we apply the almost diagonalization result to the propagators of 
Cauchy problems for constant coefficient hyperbolic and parabolic-type equations. 
We consider an operator 

m 

(56) P(d t , D x ) =d? + J2 a k {D x )d™- k , t G M, x G R d 

k=l 

where afc(£), 1 < k < m, are polynomials. They may be non-homogeneous, and their 
degree may be arbitrary. Here we set D x . = ^d Xj , j = 1, . . . , d. 
We are interested in the forward Cauchy problem 

(p(d t ,D x )u = 0, (t,x)£R + xR d 
1 d^u(0, x) = Uk{x), < k < m — 1, 

where uj~ G S(M. d ), < k < m — 1. Hence we suppose that the forward Hadamard- 
Petrowsky condition is satisfied: There exists a constant C > such that 

(57) (r, f)eCx M d , P(ir, = Im r > -C. 

This is a sufficient and necessary condition for the above Cauchy problem with Schwartz 
data to be well posed |3 1 1. Section 3.10]. By taking the Fourier transform with respect 
to x and using classical results for the fundamental solution to ordinary differential 
operators [331 PP- 126-127], one sees that the solution is then given by 

m— 1 m—k—1 



u(t,x) = ^2 dt E (ti •) * (um-i-k + ^2 aj( D x)u m -k-i-jy 

k=0 j=l 

Here E(t,x) = ^L yx o'(t,^), where a(t, £) is the unique solution to 

m 



k=l 



supported in [0, +oo) x R d . The distribution E(t, x) is therefore the fundamental 
solution of P supported in [0, +oo) x R d . 
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The study of the Cauchy problem is therefore reduced to that of the Fourier multi- 
plier 

(58) a w (t, D x ) = a(t, D x )f = J^ 1 ^, -)Ff = E(t, •) * /. 

(For Fourier multipliers the Weyl and Kohn-Nirenberg quantizations give the same 
operator) . 

Example 5.1. Here are some classical operators and the symbols o~(i,£) of the corre- 
sponding propagators for t > (a(t,£) = for t < 0): 
Wave operator — A; <r(i,£) = sin ffig| • 

Klein-Gordon operator d? — A + m 2 , with m > 0; o~(t, £) = sin ^v /47r ^ + ™ I m 

Y 47T 2 |^| 2 +TO 2 

Heat operator dt — A; o~(i,£) = e~ 471 " 2 l€l 2 * . 

We now introduce an assumption under which the multiplier a(t, D x ) falls in the 
class of pseudodifferential operators considered above. 

Assume that there are constants C > 0, r > 1 such that 

(59) (r,()eCxC d , P(ir,C) = 0^Imr> -C(l + |ImC|) r . 

It is clear that this condition is stronger than the forward Hadamard-Petrowsky con- 
dition. 

Theorem 5.2. Assume P satisfies ( |59| ) for some C > 0, r > 1. Then the symbol 
of the corresponding propagator a(t, D x ) in (58) satisfies the following estimates: 

(60) |3|V(i,0l < C (m)H+ '(a!f, (£l d , t>0, «£N d , 
TOt/t fj, = 1 — 1/r, /or a neu> constant C > 0. 

Proof. It is well known (see e.g. Proposition 1.3.2, Lemma 1.3.3], (3 1 (. Section 3.10]) 
that a(t, £) extends to an entire analytic function in the second variable, and that the 
estimates ( [59] ) imply the bound 

(61) |o"(t,C)l < e Ci(1+|ImCir , CeC d . 

Now, given £ G M d , we consider the polydisk B(£,R) = ]lj=i -R) = {( £ C d : 

10 - 01 < ^, 1 < i < d}, with i? = (1 + |a|) 1/r , a e N d . Observe that Q implies 

(62) sup \a{t,C)\<e ct ^ dR y. 
The Cauchy's Generalized Integral Formula 

d^o-(t, £ = tt^w / • • • / 77 c 77 , x a . +1 dCi • • • d(d 

(2m) d J JdBi(ti,R)x-xdB d (t d ,R) Ci - 6 Ql+ '••(&- Sd d+ 



and the estimate ( 62 ) yield 

j p Ct(l+VdR) 



(63) |^(t,OI < — < q 



ale v ' , „t(\a\+i) od 



RH ~ 1 (l + |a|) |a|/r 
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Using Stirling formula and 1/r 



1 - fi, 
1 



(l + laDMA- 



< 



a 



\a\ 



a i 



which combined with (63) provides the desired majorization (60). □ 

Observe that the assumption r > 1 in the above theorem implies < fi < 1. 
Combining Theorem 5.2 and Theorem |4.2| we obtain at once our main application. 



Theorem 5.3. Assume P satisfies ( |59| ) for some C > 0, r > 1, and set s 
max{l/2,l - 1/r}. If g £ S s s (R d ) then a{t,D x ) in ([58]) sate/zes 

(64) | (a(t, D x )ir(z)g, n(w)g) \ < Ce~^ w ~ z \ % Vz,w G 



/or some e > and /or a new constant C > 0. T/ie constants e and C are uniform 
when t lies in bounded subsets of [0, +oo). 

Notice that in the above theorem 1/2 < s < 1, so that we always obtain super- 
exponential decay, i.e. 1/s > 1 in (64). 

We now show that, if P(dt,D x ) is any hyperbolic operator, then (59) is satisfied 
with r = 1, and hence the above theorem applies with Gaussian decay (s = 1/2 in 

((64])), for windows g G Sy^(R d ). 

We recall that the operator P(dt,D x ) is called hyperbolic with respect to t if the 



direction N = (1,0, ... ,0) G 



is non characteristic for P, (i.e. its principal 



symbol - the higher order homogeneous part in the symbol - does not vanish at N) 
and P satisfies the forward Hadamard-Petrowsky condition (57). It follows then that 
the operators aj(D x ) in (56) must have degree < j and P has order m. 

The wave and Klein-Gordon operators are of course the most important examples 
of hyperbolic operators. We emphasize, however, that P is not required to be strictly 
hyperbolic, namely the roots of the principal symbol are allowed to coincide. 

For example the operator P = df — Ylj k=i a j,kd Xj d Xk , is hyperbolic if the matrix 
a^fc is real, symmetric and positive semi-definite. 

Proposition 5.4. Assume P(dt, D x ) is hyperbolic with respect to t. Then the condition 
(|59|) is satisfied with r = 1 for some C > 0, and hence 



(65) 



\(a(t,D x )n(z)g,7T(w)g)\ < Ce-^~ z \ 2 



V z,w G 



n2(/ 



if g G Sy^(M. d ), for some e > and for a new constant C > 0. 

Proof. Denote by P m the principal symbol of P and T(P,N) the component of N in 



{6e 



pd . 



■ Pm(0) 7^ 0}. It follows from the hyperbolicity assumption (see e.g. [27, 



Proposition 12.4.4]) that the symbol Q(r, £) = P(ir, ^) of P(dt,D x ) satisfies 

Q(E + tN + o-9) / i/Sg R d+ \ Imr < -C, Ima < 0, 9 G T(P,N), 



for the same constant C which appears in (57). We then deduce that 

(r, C) G C x C d , Q(t, C) = P(iT, = (Im r, Im () -CN - T(P, N) , 
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Indeed, if (r, C) = E + i(-CN-6), with 6 G T(P,N), E G R d+L , thenO-eN G T(P, JV) 
is e is small enough because N) is open, and then 

P(t, C) = P(3 - i(C + e)iV - i(0 - eN)) / 0. 

Hence — (C + Imr, Im() T(P,N), and since the cone T(P,N) is open, this implies 

-(C + Imr,ImC) -AT _ -(C + Imr) 
|(C + Imr, ImC)| ~ |(C + Imr, ImC)| ~ 
for some constant Ci < 1, which gives 

Imr > -C- 1 |ImC|- 
1 — Ci 

Example 5.5. Consider the wave operator P = <9f — A in K X M d , hence <r(£,£) = 
^~2^fff^ . Using some further information on the fundamental solution we can give a 
very precise result on the matrix decay in dimension d < 3. Namely, for the window 
g(x) = 2 d / 4 e"^N 2 (g G S^W 1 ), \\g\\ L 2 = 1), we claim that 

\(*(t,D x )MzT x g,Mz,T x ,g)\ < t e "f r-?l 2 +(l-'-l-^], e E d , d < 3, 



where (•)+ denote positive part. This agrees with Proposition 5.4 with the constants 
made explicit. 

In fact, in dimension d < 3 we know that 

a(t, D x )f(x) = j f(x- y) d/j, t {y) 

where fit is a positive Borel measure, supported in the ball -B(0, t), with total mass = t 
(see e.g. [HJ Chapter 4]). 

Hence, using T y M^ = e'^^M^Ty, 

(a(t,D x )M ( T x g,M e T x ,g) = J e' 2 ^ (M^T x+y g , M^,T x ,g) dfi t (y). 

An explicit computation shows that 

\(Tr(z)g,ir(w)g)\=e-^ w - z \ 2 , z,w G M M , 

so that 

\(a(t,D x )MsT x g,Mz,T x ,g)\ < J e~f n£'-«[ a +l*'-*-«l 2 ] d ^ y) 



which gives the claim. 

Consider now the Gabor frame G(g,A), with g(x) = 2 d / i e~' K ^\ A = Z d x (l/2)Z d 
([21, Theorem 7.5.3]), and the Gabor matrix 

Pml ,n' ,m,n — 

(a(t, D x )M n T m g, M n iT m ig), (m,n), (m,n)GA. 

We therefore have 

|r m ,y )m)n |<r OT> , W) „:=ie-f[l"'--l 2+ (l w '- m l-*) 2 + ] J (m,n), (m',n') G A, d < 3. 
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Figure [I] shows the magnitude of the entries, rearranged in decreasing order, of a generic 
column, e.g. T m / )n / j0j o (obtained for m = n = 0), at time t = 0.75, in dimension d = 2. 
In fact, the same figure applies to all columns, for T m i jn i^ n = Tm'-m,n'-n,ofl- This 
figure should be compared with [6, Figure 15], where a similar investigation was carried 
out for the curvelet matrix of the wave propagator on the unit square {d = 2) with 
periodic boundary conditions. It turns out that the Gabor decay is even better, in 
spite of the fact that we consider here the wave operator in the whole R 2 . 




10 10 10 

Coefficient number 



10° 



Figure 1 . Decay of a generic column of the Gabor matrix for ^2prDT^ 

in dimension d = 2 and at time t = 0.75, with window g(x) = \/2e~ n ^ 2 
and lattice 1? x (1/2)Z 2 . 



We now present a class of examples of operators which satisfy Theorem 5.3 and are 
not hyperbolic. 

Example 5.6. Consider the operator 

(66) P(d t ,D x ) = d t + (-A) k , 

with k > 1 integer. When k = 1 we have the heat operator. Its symbol is the 
polynomial 



We claim that it satisfies (59) with r = 2k. 
If r G C, £, rj G R d and P(ir, £ + irj) = then 



Imr = (27r) 2k Re(\e + 2i£-r,- \ V \ 2 ) k = (2^1^ + Q(£, v), 



2\k 



\2k\t\2k 
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where Q(£,r]) is a homogeneous polynomial of degree 2k, with 0) = 0. Hence, for 
every e > and some constant C e > 0, 

\Q(t,v)\< E \Z\ j \v\ l <^\ 2k + c e \ v \ 2k , 

j+l=2k,j>0,l>l 

where we applied the inequality ab < + ^ b ^ q (with a = |£p, b = \r]\ l , p = 2k/ j 
and q = 2k/l) to the terms of the sum with j > 1. Taking e = (27r) 2fc /2 we get 

(2vr) 2fc 



Imr > 



\e k -c\r,r>-c\ v r, 



\2k 



2k 



which proves the claim. 

Hence, for the operator P in (66), Theorem 5.3 applies with r = 2k and s 
1 - l/(2fc). 



Remark 5.7. Theorems 4.2 and 4.4 consent also some applications to equations with 
variable coefficients. As elementary example consider the transport equation 



dtu - i (Ej=i a 3 D *j + E?=i u = 
«(0,x) = uo(z), 



with a,-, 6,- G M. The fundamental solution of the problem has symbol 



a 

a(t, x,£) = exp \it ^^(ojCj + bjXj) 



satisfying (44) with m = 1, s = 0, and the preceding arguments apply. More generally, 
it seems possible to treat similarly the pseudodifferential problems of the form 

dt — ia w (x, D)u = 0, u(0, x) = uq(x), 

where a(x,£) is real- valued and d x ^a(x,^), as well as higher order derivatives, are 
bounded in M 2d , satisfying suitable estimates. 
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